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Abstract 

We prove the invariance of the contact angle in liquid-solid wetting 
phenomena : an electrified droplet is spreading on a solid surface. The 
drop is minimizing its energy. We express the differential of this energy 
with respect to the shape of the drop and deduce necessary conditions for 
optimality . By variational method, using well-chosen test functions, we 
obtain the main result about the contact angle between the drop and the 
solid. 
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Introduction 

In this article we give a proof of the invariance of the contact angle in a solid- 
liquid wetting phenomenon. To model this experience, described below, we 
use an energy minimization with respect to the shape and we obtain the main 
result by variational techniques. The phenomenon we are dealing with is named 
electrowetting. 

Electrowetting is a technique allowing to modify the affinity between a 
solid and a liquid, by the introduction of an electric field. This phenomenon 
has been discovered in the XIX*^ century by Gabriel Lippmann on a mer- 
cury/electrolyte solution system ([ISj). It has been then studied on more general 
systems ([HIIS]). We can find today many industrial appHcations. Particulary 
in optic with variable focal lenses (see the web site of the society of B.Berge : 
http://www.varioptic.com/en/ and [2j), and pixels for electronic paper (^9.i20j). 
It is also use in microfiuidics, for example in chip design (|18]) and has biomed- 
ical applications (jllj) 

Consider a liquid drop on a polymer film. The equilibrium of the drop results 
on superficial tension forces and gravitational force. The drop shape is a quasi- 
spherical calotte with a contact angle between water and solid which is given 
by Young's angle (1805)^. We applied then a constant voltage 0o between this 
drop and an electrode placed under the insulating polymer. The charged drop 
and the electrode create a capacitor. A simple model considering the system 



1 



as a plane capacitor predicts the total spreading of the drop when the voltage 
becomes higher. However physical experiences show a locking of the contact 
angle when the drop is bound to a voltage higher than a critical value. 

Many explanations have been proposed for this saturation of the angle. Most 
of them give a crucial role to the divergence of the electric field in the vicinity of 
the triple line (solid- liquid-gas interface). The understanding of this phenomenon 
is slowed down by the misunderstanding of the geometry of the drop near the 
triple line (or wetting line). 

We prove in this paper that the contact angle is independent of the applied 
potential and that the value equals Young's angle, as observed in [5]. 

The problem will be described in the first part and the equations of the 
model established. 

In the second part we will recall some shape optimization results adapted to 
the model. 

In the third part we will establish necessary conditions for optimality which 
will be exploited in a fourth part during the calculation of the contact angle 
value. 

1 Description of the problem 

We study the evolution of the shape of a droplet located on a substrate and 
subjected to a constant voltage. A 3D model (as found in |4j) allows us to 
examine axisymmetric case and envisage developments to non symmetric shape 
for high voltages. 

1.1 Hypotheses 

Assumptions : 

i) The applied electrical potential (fio is continuous 

ii) The liquid drop is a perfect conductor. 

iii) Electrostatics effects are negligible far away from the drop. 



1.2 Experimental device 

We consider an orthonormal basis in M'^. The top side of the polymer film, 
where the droplet is posed, is the plane [Oxy). We use the indices L,S and G 
to refer to liquid, soHd and gas domains. A couple of indices LS,LG,... relates 
to a liquid-solid, liquid-gas... interaction 

Vie is the bounded domain where the experimental device takes place and 
where calculations are directed. 
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We denote by : 

flL the liquid domain and Tl its boundary. 
Tls the hquid-sohd interface. 
Tlg the Hquid-gas interface. 

the gas domain and Tq its boundary. 
fls the sohd domain and Fg its boundary. 
To the boundary of fig where the counter electrode is applied. 
FeG and T^s the other external boundaries of fig 

EG, £S and El permittivities of ilc, and Ql, respectively. 

Ni, i equals to L, G, S, LS, LG, GS..., normals to the surfaces Fj, respectively. 



1.3 The electrostatic model 

The apphcation of an electrical potential (po between the counter electrode Fq 
and the drop flL creates an electrical potential in the entire space. The drop 
is supposed to be perfectly conductive and the potential is also constant in Ql- 
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But the charges distribution is not constant : it depends on the shape of the 
drop i.e. on fig \ f^L- 

(fp' is the solution of a system of partial differential equations : 

: G, S. 



div(eiV0f ) = 





in 


</'g = '? 


^0 


on Tlg 


4>% = ^ 


^0 


on Tls 


^^s = 





on To 



At the solid-gas interface, wc have the following transmission relations : 

^ = 0g ^ on TsG 
EG V(/.g .TVg = -es V<^g .iVs on Tsg 

On the artificial boundary, we impose : 

eiV(t>f.Ni=0 onTei, i = G,S. 

This system which gives the potential can be rewritten in a weaker form. 
We set : 

HoiCl) = {V e H^{n); = on To U Tls U Tlg} 

H{Q) = {V- G H\n);^ = on To, = (Ao on Tls U Tic} 
Denoting £ the map defined on f2 by : 



Eg on r^G 
£s on ^5 



The weak formulation is : 

find (jP e such that 



iFV)n 



where < ., . > is the scalar product of M^. 
In the following we denote : 

an{(l),tp)= e<V(l),Wtp>dQ (1) 
Jq 

Thanks to Lax-Milgram's theorem we are able to prove that this problem 
admits a unique solution (fP G H{Q). 
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1.4 The energy of the drop 



The potential (/>o > being given, the equilibrium shape of the drop flL cor- 
respond to a minimum of the energy of the system. We take into account the 
superficial tension force, gravitational force and electrostatic force. We denote 
o'LS, o'GS and (Tlg surface tensions relative to the different interfaces. 
For a drop fii and a potential 0o, the energy of the system is : 

6{^L,4>o)= P9 / zdfl + / aLoda + / aLsda + / acsda 

•f^L '^^LG '^^LS '^^GS 
V ' ^ V ' 

Potential energy Capillary energy 

e\\/(j)^^\^dn 
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Electrostatic energy 
Up to an additive constant, we have : 



£{p.LAo)=P9 zdVL+ f aLGda+ I {(JLS-crGs)da-]- e\W(tP\^dn 

where the electrostatic energy is affected by a minus sign because this energy is 
imposed by an exterior generator. This integral depends on fi^ since O = Qe\^L 
and 4>^^ is the solution of the variational formulation {FV)a. 



1.5 Search of the optimal shape 

The drop of volume V is submitted to a voltage (po- The optimal shape fi^ 
of the drop correspond to a minimum of the energy : 

£{^ll,(j)o)= min £{^1l,(Po) 

{nj:,;Volume(ai,)=y} 

The evaluation of the function £{Ql,4'o) requires the resolution of a partial 
differential equation's problem on the domain = ile\^L- ^^s is fixed but i^c 
and also O depends on : It is a major difSculty of the problem. 

Obviously it is equivalent to give Oi, or to give O, and in the following we 
take O as a variable. 

By introducing the parameters : 

pg CTLS - CTGS c 1 

a= , /u= , 6 = 



ctlg ctlg ctlg 

our problem is equivalent to minimizing 

J{n)=-a I zdn + ii I da+ I da - ^ [ e\V4P\^dn 
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We want to minimize J on a set of admissible domains obtained by small 
smooth deformations of a reference domain. This point will be clarified on 
paragraph 2. 

The goal is then to determine a necessary condition for optimality for a 
domain fi. That's why we are going to give a sense to the derivative of the 
functional J. The main difficulty is that we want to derive on a set of domains 
which has not the usual stucture required to define a derivative in the ordinary 
sense. 

2 Necessary conditions for optimality 

To obtain a necessary condition for optimality, let us make precise the class of 
domain on which we are minimizing. 

The theory we use is detailed in [16], [lO] and [21] . 

2.1 Admissible domains 

We need to give a sense to integral formulation on domain or on boundary of 
domains. We will work with open sets fl with Lipschitz boundary. For this 
reason, we choose to take deformations of open sets with Lipschitz boundary, 
obtained by sufficient smooth maps, in order to keep the lipschitzian feature of 
domains ([21]). 

• We denote 

with the infinity or sup norm 

||C/||oo = sup|C/(a;)| +sup|£iC/(a;)l 

where DU is the differential of U 

As we only consider bounded domains, for all U in C^(ri,M'^), we have U 
and DU uniformly bounded on Cl. 

• We introduce the set of admissible displacements : 

UiW^) ^{Ue Ci(n,K3); \\U\\oo < 1, U.^ns = 0, and C//r, = 0} 

Finally for U G Ci(n,R3), n + U denote the set {Id + U){n) 

• Let be a fixed reference domain of the type described in paragraph 1. 
Define 
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We are searching a necessary condition for optimality for the solution of 
the following problem : 



where C{Q) is the volume constraint, more precisely : 
C{n) = \o\{ne \n)-V, recalling that = \ 

We are now going to give a sense to the differentiation of the function J 
comparatively to a domain fl with the concept of shape derivative. It is a 
classical notion which can be find in detail in [16j. Here a weaker version 
is given, which is still sufficient for our problem. 

• Directional derivative. 

Definition 2.1 A function J defined on Vad and with values in K. has a 
directional derivative at a point CI of Vad in a direction U G if the 

function J*: V ^ J{{I + V){il.)) defined on U {CI, M.^) with values in 
M has a directional derivative at the point in the direction U (in the 
usual sense in C^(M'^, M"^)/ The directional derivative of J at fl in the 
direction U is denoted : DJ{n).U := DJ*{0).U. 

We are now able to write a necessary condition for optimality for J. 

For fl € Vad, and A G E, we denote A) = J{fl) — XC{fl) the lagrangian 
of our problem of optimisation under constraint. We are going to find a saddle 
point of £ O)- We are searching a necessary condition for optimality for a 
couple (f2*, A*) to be a saddle point of £ with fi* e Vad and A £ K. 

2.2 Necessary condition for optimality 

Proposition 2.2 // is a saddle point of C, then if J and C admit 

a directional derivative at fl* in the direction U e U{n* ,M!^) , DJ{fl*).U = 
X*DC{Q*).U. 

Proof : {n*, A*) saddle point of £ e Vad, VA £ R, 

c{n*,\) < c{n*,x*) < c{n,x*) 

Let U G U{TF,M.^). For t e [-1, 1], tU is an element ofU(n*,M.^). And so the 
set ft* + tU is an element of Vad too. 
Then we have : 




£{n*,x*) < c{n* + tu,x*), vte[-i,i] 
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By making explicit the values of £, we get : 

j{n*) + x*c{n*) - j{n* + tu) - \*c{n* + tu) < o, vt e [-1, i] 

By taking definition's notation, we deduce that : 

J*{tU) - J*(0) + A* [C*{tU) - C*(0)] > 0, Vt e [-1, 1] 
Taking i > 0, we obtain : 

t t - > L > J 

Let t tend to 0, we obtain : 

DJ*{Q).U + X*DC*{Q).U > 

With the same argument but with i < 0, we finally obtain : 

yU eUiTF,R^), DJ*{0).U + X*DC*{0).U = 
That is to say, by definition of shape derivative : 

vc/ e w(IF, ]r3), Dj{n*).u + x*DC{n*).u = o 

■ 

We are now searching the expressions of the derivative of J and C with 
respect to a domain. 

3 Derivative of the drop's energy with respect to 
its shape 



Let ^l* C "Dad, such that 



J(n*) = min j(n) 

C(Q) = 



where 



J{Q) = -a [ zdn+ [ da + ^i [ da-^ [ e|V(/>"pdl7 



ivt^ ^ ^ ^ 

Jgrav Ji^o Jls Jdi 

and 

c{n) = Voi(r2e \Ti) ~v 

J is the sum of four terms Jgrav, Jlg, Jls, Jei- The first three terms of J 
and C are integrals on a surface or a domain of a fonction independant of this 
domain. We have classical results on the shape derivation of such terms, (see 
for exemple jlGj). So we won't give more details for the derivation of this terms. 

We are going to spend more time on the singular term of our problem : the 
electrostatic contribution Jei- 
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3.1 Derivative of gravitational, capillary and volume con- 
straint term 

Let n £ Vad and U € U{Ti, R^). 

The derivatives at the point ft and in the direction U is obtained by fohowing 
Definition 12.11 and the results of p2j : 

DJgravii^)-U = Dj;,^^{0).U = a [ U,dn + a [ zdiv{U)dn (2) 

• In the same way, we express the differential of C which looks like the term 

J grav • 

DC{n).U ^ [ diy{U)dn (3) 

• The derivative of the terms Jlg and Jls are obtained hkewise using a 
result about derivative with respect to a surface : 

DJLom.U = / div(f/)da - / < N^,' DUN^ > da (4) 

DJLs{^).U = iil div(f/)dt7 - / <N^s' DUN^s> da (5) 
where *DU is the transposition of the Jacobian of U 

3.2 Derivative of the electrostatic energy 

We are studying much more in detail the electrostatic term Jei ([15]). 
Let V e U{Ti,R^), by Definition ([III]) 

^ Jn+v 

To lighten the notation, let us pose F :— Id + V. 
As y e U(Ti,R^), F is invertible. 

• G H{F{Q j) is the solution of the variational problem {FV)p(^Qy 

• The map Tq : ip £ Ho{F{fl)) t-^ ip o F £ Ho{i^) is an isomorphism from 
Ho{F{n)) to Hoin). Likewise we define T from H{F{n)) to H{n), because 
00 is a constant ; T is an isomorphism too. 

We consider the two transported variational problems : 
r Find £ H{F{n)) such that 
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(FV) 



Find (j)^ e H{n) such that 

y-ijj e Ho{n), e <* {DF~^ o F)V(j)^,^ {DF'^ o F)Vil) > | dei{DF)\dVL = 



Problems {FV)F{n) and (-F^^)f2 are equivalent. 

In the following we note, for 4> G H{VL) and -0 £ i?o(f^) 



a{F,(P,ij) := [ e<* {DF'^oF)\/(l),* {DF-^oF)^!!; > \det{DF)\dn (6) 

Jf2 



By unicity of solutions of the two variational problems 
and so 









/f(o) 




Jn 



Thus, let us pose for G H{Q), 

jli{F,^)^^ f e\'iDF-^ o F)^(t>\^ \det{DF)\dn 



(7) 



In the expression of J*;, variables F and (j)^ have been decoupled. 

By the implicit function theorem we can show that Jei is and that 
V ^ is differentiable ([21]). 

A classical optimal control result ([l6j V-10) allows us to write 

yU G U{n,R^), DJli{Q).U = DFJei{Id,(f>^^).U - DpHld, (tP,p^).U 

where Dp denotes the partial differential with respect to the variable F . 
G Ho{Q) is the adjoint state solution of the equation 

V0Gi?o(f^), D^d{Id,cb^,p^).(j)^DMld,(b'^)4 (8) 
The differential of ((6]) and ([7]) with respect to (p when F = Id gives : 

y(b£Ho{n), Dfa{Q,ctP,p'')4^an{p'',(t>) (9) 
and (jP" being the solution of the variational problem (FF)n, 

D^Jeiild, = ~5an{(iP, 0) = (10) 

We deduce of dH]), I© and (HO]) that V0 G Ho{n), an{p^,4>) = 0. 
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And so as 6 Ho{VL), 

From this we deduce, 

VC/ eU{n,K^), DJ*i{Q).U = DFJli{Id,<j>^^).U 

Now we look at an integral on a fixed domain and using the differentiation 
under the integral sign formulas (|21|). we obtain finally 

DJ^i{n).U = DJ*i{0).U = [ e|V0"pdiv(t/)dn 

S f 

+ - e <{^DU + DU)V(tP,V4P >dn (11) 
2 Jo. 

In the following, we simplify the notation by denoting Lgrav, Leant, L^s, 
Llg and L^i respectively for DJgrav, DC, DJls, DJlg et DJ^i. 

3.3 Formulation of the necessary condition for opti- 
mality 

By proposition (|2.2p . if a pair (O*, A*) of Dad x M is a saddle point of C 
then 

Lgrav{fl*).U + Llg{^*)-U + Lls{^*).U + Lel{n*).U = \* Lcont{n*) .U 

(12) 

Terms used below are defined by (g]), (HJ and (fTTI) . 

The formulation obtained there is verified for all 3D optimal domain ofVad- 
It can be used for numerical simulations particulary by high potential (/jq. 

We are giving a first application to calculate the contact angle for axisym- 
metric shape. 

4 Calculation of the contact angle for an axisym- 
metric shape 

This paragraph contains the proof of the main result : in the case of an axisym- 
metric optimal shape, the contact angle is the static Young's angle. 

First, we will write necessary conditions for optimality introducing the ax- 
isymmetry in the model. 

Then, we will by a judicious choice of direction of deformation calculate the 
value of the contact angle. 



= 
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4.1 The axisymmetric problem 

Let's suppose that the domain f2 of K'' is axisymmetric. We choose to express 
a point of the space in cyhndric coordinates. 

In an orthonormal basis of we denote ui the 2D domain associated to O 
in 3D. We define in a similar way Wl, ws, wg, Wg, 7ls, Jlg, IGS, 7e- 
All are defined as in the figure : 



leG 




7o 



7e is constituted of external Hquid (7eL), solid (7es) and gas (7eG) boundaries. 
A is the contact point, and Qa the contact angle. 

4.1.1 Notations 

For w , we define the analogous spaces to those defined for the domain O. 

• U{u!,R'^) = {ue Ci(w,R2); ll^ll^ < = 0, et u/^^ = O} 

Let us suppose that we are at the minimum of the energy fl* and consider 
Lo* the associated 2D domain. 

• We choose a direction of deformation U € U{ft*,M.^) which is invariant by 
a rotation around the z axis. 
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We are able to find 



(r, {ur{r,z),Uz{r,z)) 

such that u G ll{uj^, R"^) and if 

U : IF ^ R3 

(r, 9, z) 1-^ {ur{r, z) cos{0), Ur{r, z) sin(0), Uz{r, z)) 



then 

U{rcos{e),rs.m{e),z) = U{r,e,z) 

• In the same way, we have analogous relations and notation for the different 
normals to the boundaries. When there is no ambiguity on the considered 
surface 7 (resp. F), we will denote it (resp A'^) for (resp Np). 
It exists 



{r,z)h^ {nr{r,z),nz{r,z)) 

such that if 

{r,6,z) {rir (r, cos(^) , (r, sin(^) , Hz (r, 2;) ) 



then 

N{rcos{e),rsm{e),z) =N{r,e,z) 

The deformation direction and the normal are also entirely determined by 
the maps of the plane u and "n . 

• We have 

div(?7) = y + div(u) (13) 

and 

< N* DUN >=< 7?,* Du'n > (14) 
for the normal to the different considered surfaces. 
We will note 

div-y(u) = div(u)— < li* Du'n > 
the surface divergence of the surface 7. 
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We define a 2D potential associated to the 37? potential too. 
We will note 

Muj*) = e Hl{uj*);u = sur 7o* U ^ls U iIg\ 

Hiuj*) = {u e Hl{uj*);u = sur 73, m = 0o sur U 72^} 

where Hl{uj*) is the Sobolev spaces for the measure dv = rdrdz. 

We will note the analogous spaces for the case with the measure v. 

Consider the following variational problem 

Find ip"^' e H(uj*) such that, 



(/«) 



WveHoiuj*), J^, e <Vtp'^' ,Vv > rdrdz = 



where < ., . > is the scalar product of R^. 

In the following we will denote for ip e H{uj*) and v e Ho{uj*), 

a^»{ip,v)— / er < V(^, Vw > da; (15) 

To <jP is also associated (/^'^ the solution of the 2D problem given by his 
weak form. 

We can verify that (fP' is a weak solution of the problem which gives the 
potential on Q.* if and only if Lp'^ is a weak solution of the associated 2D 
axisymmetric problem. 

• Let us precise the boundary ^lg in the vicinity of the triple point. By 
definition of admissible domains, w* is a Hpschitzian open set. Let A be 
the triple point. If we suppose that the contact angle at A is in the interval 
]0, 7r[, then it exists an open ball Ba centered on A of radius do and a 
function / : M ^ M such that, (r, z) G 72^ n Ba <^ r = f(z) (Implicit 
function theorem) . With such a parametrization we can express the cosine 
of the wetting angle at ^ = (/(0)i 0) '■ 

cos{Qa) = , ^'^^^ (16) 

4.1.2 Formulation of the necessary condition for optimality in the 
axisymmetric case 

By adopting the previous notations, and by noting Igrav, Ilg, Ils, hap and Icont 
the associated terms in 2D, identity lfT2|) becomes 



3A* G M tel que Vw e U{uj*,R^), 

1{UJ*).U := lgrav{0J*)-U + Ilg{^^*)-U + lLsi^^*)-U + lel{^*)-U = \* IcontiuJ*) -U 

(17) 

We express those differents quantities in detail in the following paragraph. 
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4.2 A particular choice of admissible directions 

Suppose that we are at a minimum of energy w* 
We know that at this minimum 

3X* E R such that Vu e UiZJ^,^^), l{uj*).u = X*lcont{uJ*).u 

We want to use this equality in order to extract an information on the 
geometry of the drop in the vicinity of the triple point A. The idea is to find 
particular deformation's directions u^, p G N* whose support, centered at A 
focuses on A as p tends to +00. 

The sequence (M^)pgN' which is chosen is defined by : 

uP : M2 ^ 

{r,z) ^ K,0) 

where 
uP.: R2 ]r 

f ^ J exp(^y^ I -) if (r-/(0))2 + z2 < 1 

else. 

The support of uP is the ball of center A and of radius - . It is included in 

P 

the neighbourhood Ba where j^q is parametrated by /, if p is higher enough. 
We have e M^). 

In the following paragraph, we study the limit ofl{Lu*).uP and of lcont{^*)-uP 
as p tends to +00. 

4.3 The limit of 

In the following, we denote Bp the open ball of center A and of radius - and 

P 

dBp its boundary. 

Let us study each terms which appeared in l{uj*) 

4.3.1 The term lgrav{uJ*)-uP 
With ^ and ([13]), we have 

lgrav{i^*)-u^ — OL I uPrdrdz + a / zuPdrdz + a / zdiw{uP)rdrdz 

It is clear that 

/ ruPJx = 0, since < = (18) 
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As \u^\ < 1, we have 



/ zUrdrdz < j \z\drdz < n— 



(19) 



duP, 



As uf = 0, we have / zdiv{uP)rdrdz = / z—-^{r,z)rdrdz 



By posing R := ^^"^TWF+i^ 



1(0) 



and the upper bound 



dr 



(r, z) 



{p^R^ ^ 1)2 
R 



exp 



1 



< 



1 



R 



From the study oi g : R > 
C 

\g{R)\ < ~j where C is a positive real, 
and so 



exp 



(p2i?2 „ 1)2 ^ yp2^2 _ 1 



|div(u'')| < 2Cp 



We finally obtain 

z(]xw{vF)rdrdz 
Then with ifTO]) and IpTj) . 



< 2Cp 


/ zrdrdz 







we deduce that 



(20) 



= 2C^(/(0) + i) (21) 



lim lgrav{i^*)-uP = 

p — ^+oo 



(22) 



4.3.2 The term Ilg{^*)-uP 
From ig]), (dHl) and ([H]), we deduce 



Zlg(w*)-w^= / uP{s)ds+ / div-y(uf)rds 



Let 7(p) be the intersection ordinate of Cp with 7|,(3. As p > po, we use the 
parametrization of "f^Q by /. We have 



uP{s)ds 



7l,G 



7(p) 



so 







/ uP{s)ds 


-I 







uUf{z),z)^i + r\z)dz 

^1 + r\z)dz = l{^la^Bp) 
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and then 



lim / uP(s)ds = 0. 



(23) 



The study of the term div-y (u*')rds which contains the surface diver- 
gence is a httle bit more fastidious and contributes to the final calculation of 
the contact angle. 



We show that the hmit of this term is 
to +00. 



/(0)/'(0)(l + /'2(0))-| 



as p tends 



/ div-y(u^)rrfs = / div(u^)rds — / <* DvFli; 'n > rds. 
•'^LG •'iIg •'iIg 



"Ilg 

As 



and n 



V dz 



i 



Therefore 
/ div-y(M^) rds 



dz 'l + /'2(z) 

n(p) , rhiP BiiP 

f{z)f\z){l + f\z))-Hnzf-^ + ^]dz 

nip) , 

f{z)nz)ii+r\z)r-^g-^m{z),z)]d. 



fizWl + f^z)dz 



Let us pose hiz) := + f'^{z))-i. 

We are now showing that 



lim 

p— ^+00 



div^K)rds = — h{0) 



As / belongs to C^([0, 7(p)], M), h is a continuous function and so particulary 
continuous in : 

V£ > 0, 35{e) > such that if \z\ < d{e), then \h{z) - h{0)\ < e 

We have : 



1^ h{z)-miz),z)]dz + -h{0) 



'lip) 

< I \hiz)-hm 
'0 



dz 
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pip) Q I 

because / —[u^if iz) , z)]dz — — 
Jo oz e 

As 7(p) tends to as p tends to infinity, 

Ve > 03p*(e) e N* such that Vp > p*,7(p) < (5(e) 



By the same reasoning as the one used to establish l|20p and as /' is bounded, 

az 



we are able to show that 



Moreover, as 7(p) < -, we deduce that 
P 



< Dp, with D a real constant. 



nip) f) 1 

hiz) — [uP{f{z),z)]dz + -hiO) 



< e 



7(p) 



g-J<if{z),z)] 



1 



< £Dp- < De 
P 



dz 



To summarize 

Ve > 0, 3p*{e) e N* such that Vp > p,(e), 

d 
oz 



(M^)-K(/(z),z)] 



-h{0))dz 
e 



< e 



We have also proved that 



lim 



rdiv-^(uP)d.s = --h{Q) 



(24) 



4.3.3 The term Ils{uj*)-vP 
By (HI), (Uni and (HI]), we have lLs{i^*)-u^' ^ fJ^ u^ds + /i / diVj{uP)rdr 
As 7LS n support(M^) 



/(o)--,/(o) 
p 



and is bounded by 1, we have 



<ds < (/(O) - (/(O) - ^)) = ^ 



so 



lim 

p — >+oo 



Ids = 



(25) 



7iS 
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It is the second term containing the surface divergence / divj{u^)rdr on 



We have *DuP 



7^5 which is going to contribute to the calculation of the contact angle, as for 
the term Ilg- To estimate this term, we have to express *Z?uP and ~n . 
( ^< 

and ^ = I 



I Ahi^lvF)rdr — I -—^{r,0)rdr 

J lis JfW-j; * 



Integrating by part, we find 



div^K)rdr = [r<(r,0)];[[;;;_, 



/(o) 



fioy 



<(r,0)dr 



The second term on the right side of this identity has a null limit as p tends 
to +00 because is bounded by 1. 

As <(/(0),0) = i and <(/(0) - ^,0) 0, we have 



and thus we conclude that 

/(O) 



Hm / div^ (u^')rdr = (26) 



'7 

p — >+oo ' 



4.3.4 The term lei{uj*).uP 
From ifTTj) and ifTS)) . we deduce 



leiiu;*).uP 



4:1 suP\Vip'^'\^drdz-^ f £|V^"*|2div(Mf)rdrdz+ 

^ [ e < CDuP + DuP)Vip'^' ■ V^"' > rdrdz 

These terms are containing the electrostatic contributions in the derivative 
of the energy. 

For the first term / euP\yip'^ \'^drdz, we use the dominated convergence's 



theorem of Lebesgue. Indeed, we have 

^KxbJ'^v'^ P ^ p.p., as p +00 
Moreover, we can find a positive constant M such that 
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As we know that (^"* G Hl{uj*), we deduce that |Vi^'^*p G 

The dominated convergence's theorem of Lebesgue allows us to claim that 
lim / euP\Vip'^'\'^drdz (27) 



The second term / s| V^- Pdiv(.^).d.d. which appears in is more 

fastidious to treat. It is necessary to compare the behaviour of the two singular 
terms in the vicinity of A : div(uP) and jV^?'^ p. 

We have an upper bound of the divergence of given in (|20|) . For | V(y9'^ 
we can make precise the behaviour of its norm in the vicinity of the point A 
in each domain and loq. 

We choose to work with polar coordinates (p, 9) centered at A. 




Adapting a theorem due to K.Lemrabet ([l2]) which split the potential into a 
regular part and a singular part in the vicinity of the edge A, we obtain 

Theoreme 4.1 It exists a unique c G K. such that for i G {G, S} : 

- cS. e HUiot) 

with 

^) = n ^^ P^ - - ()A)))vip, 0), e G]0, tt-OaI 

cos(A(7r — Ba)) 

Ss{p, 6) = ^^p' sin(A(0 + 7r))r;(p, 6), G] - ^, 0[ 
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where : 

■q is an infinitely continuous function with a compact support, such that rj = 1 
in a vicinity of A and becomes zero outside of a ball centered on A. 
A is the unique solution in ]0, 1[ of the equation 

e5tan(A(7r — Qa)) — — £Gtan(A7r) 



By the equation verified by A, we have A e]-, 1[. 

We are now giving an upper bound of the Ll norm of the modulus of the 
gradient of ipf on u* n Bp denoted \\V(pf W^l^ (^uj'nBp)- setting gi = ip"^ — 
cSi e H'^{uj*), we deduce 

llV^f |li2(^.nB,) < 2 [lIcV^.llijKns,) + \\^9^\\luu.tnB, 

It remains to find an upper bound to the terms which appear in the right hand 
side of the inequality. 

We choose p big enough such that ri/B^ = 1 and we consider p> p. 

Calculating the explicit expression of VSi, we obtain the following upper 
bound : 

where D* is a constant. 

Then the upper bound of its Ll norm in the vicinity of the point A is 

l|v5dli^(..nB,) < (/(o) + -) / \ys,{p,0)\' pdpde 



<D*{f{Q) + -) f p^^-HpdO 
P Jb, 

= D*{f{Q) + -)^p-'^ 
P A 



(28) 



Using the Sobolev's injection in three dimension and in axisymmetric case 

m p-27) 

we obtain an upper bound of 1112(^^.^3 ) by Holder inequality : 

/ \^g.,\^du^( XB, \y9^?dl,<{^{Q) + l)\\xB^\^i,,^\y9^\\ll( 
•luj'nBp Jul' — V — ' P ^ - 

i?K)i?K) mes(Bp)t 
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From the two estimates ([28]) and l(29|) . we deduce : 

l|V<^r'lli^(.-ni3„) < 2(/(0) + l)[D*^p-^^ + K,mes{Bp)l_ 



for i e {S, G} 

We can also write an estimation of the gradient norm of the potential on the 
whole domain uj* 

l|V^"*|li.(^.nB,) <2(/(0) + i) [d*jp-'^ + K*mes{Bp)l] (30) 

where K* = Ks + Kg- 

Now we are coming back to the term / s\Vf'^ diy {u^)rdrdz. For every 
p > p, we have 



e\\/ip'^'\^div{uP)rdrdz < MDp / \Vip'^' \^drdz par ([201 



< 4MDp(f(0) + -) 
P 



A pn 



and as 

-2A + 1 < 0, 



we deduce that 



p 

It remains to examine the term 



lim / £|V(^'^*pdiv(M^')rdrdz = (31) 



f e< CDuP + DuP)Vip'^' ■ Vlp"^' > rdrdz 



This term is of the same order as the previous term. 

dz 





We have 'DuP + DuP ^ g^j; 9/ 



dz 



Denoting V(^'' 
we obtain 



or oz 
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from what we deduce 

e <* DuP + DuP)Vip'^' ; Vtp"' > rdrdz 



< 2M 



or az 



rdrdz 



< 2M 



l^WVip'^ \^rdrdz + - |^||V(p" iVdrdz 



•nB„ 



dr 



•nB„ 



dz 



We have 



duP 



dr 



<2Cp, 



duP 



< Dp and r < /(O) 



By an analogous reasoning, we can show that 

Urn / e <* DuP + DuP)Vip'^' ; Vip'^' > rdrdz = 



(32) 



4.3.5 Conclusion 

In the same manner, we can prove that 



Hm lcontiuJ*).uP = 

p — *+oo 

From we deduce 

hm l{u;*).uP = -[-h{0) + ^/(O)] 
p— >+oo e 

4.4 Value of the wetting angle 

At the minimum of energy uj* , we verify : 

3X* G M such that for all p big enough, ^(1^*).^^ = X*lcont{^^*)-uP . 

So 

lim 1{lj*).uP = A* lim lcont{^*)-uP 

p — ' + 00 p — >+oo 



(33) 



(34) 



From l(34|) we deduce that 
And so if /(O) ^ 0, we conclude : 



h{0) = /i/(0) 
fjO) 



But we know that the right hand side is in fact the cosine of the contact 
angle 

cos(8a) = — ^ 

By the definition of fi and Young's angle, we deduce that the value of the 
contact angle for an optimal drop with an axisymmetric geometry is Young's 
angle and this is vaHd for all values of potential appHed. 
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5 Numerical work in progress 



We recall that a simple model consists in considering that the system is a plane 
capacitor. In this case, the contact angle is given by the relation : 

cos(6i) = cos(6iy) + -^^V^ 
2eaLG 

where e is the thickness of the insulator. It predicts a total spreading of the 
drop on the polymer. In fact we can admit that it gives a value of a macroscopic 
contact angle. 

We use two numerical approaches. The first one is a macroscopic one and 
the second one a local one. 

From the code "Electrocap" developed by J.Monnier and P.Chow Wing Bom 
([TS]. ^^), we introduce a treatment of the singularity of the potential. In order 
to compute the singularity of the potential we use the Singular Complement 
Method as presented in [6j. 

This numerical approach gives same qualitative results as in pjtj (which 
is without the treatment of the singularity). It shows a deviation from the 
shape predicted by the plane capacitor approximation : the contact angle is 
higher than predicted by the plane capacitor approximation and the curvature 
increases sharply near the triple line. But the effect of the treatment of the 
singularity seems to be deleted. Furthermore it doesn't show that the contact 
angle is constant. 

We present numerical results for given volume and physical parameters. 
In Fig[T] we present the macroscopic shape of the drop obtained for different 
voltages. 



200V 
400V 
600V " 
800V 
900V 







Figure 1: Left, Macroscopic shape of the drop for a Young's angle of 90° ]Right, 
Macroscopic shape of the drop for a Young's angle of 60° 

A second numerical work is also in progress. The global approximation 
(using the treatment of the singularity) is used together with a local model 
(given by an ODE in the vicinity of the contact point). 
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In Figdlare presented values of the contact angle obtained by this numerical 
approach. Values are compared to the plane capacitor approximation. 




Figure 2: Left, Value of the contact angle for a Young angle of 90°; Right, Value 
of the contact angle for a young angle of 60° 

The numerical contact angle appears to be constant as theoretically pre- 
dicted in the last section and proposed in [5]. The numerical curvature appears 
to explode near the contact point. In Figl3]we present the curvature value for 
points of the drop at 500^ for a Young's angle of 90°. For other voltages and 
values of Young's angle, the qualitative behaviour remains the same. 



at 500V for a Young angle ot 90° 



Curvature values at 500V for a Young angle of 6( 



1 0.8 • 




Figure 3: Left, Value of the curvature for a Young angle of 90° at 500V; Right, 
Value of the curvature for a Young angle of 60° at 500V 

All this results are in accordance with the theory and experimental works 
(see e.g.[3]). 
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6 Conclusion 



We have proved that the contact angle in electrowetting remains constant for 
all potential applied and it equals Young's angle. It is the result that has been 
predicted in [5j. A numerical work is also in progress. Because of the singularity 
of the potential in the vicinity of the contact point, it uses a local model near 
the triple line. The contact angle computed appears to be constant. 
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